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1 Introduction 



Cellular Automata (CA) were first introduced in the early fifties by J. Von Neu- 
mann pQ in his investigation of " complexity" , following an inspired suggestion 
by S. Ulam. S. Ulam himself was among the first to study intensively CA as 
well the nice properties of simple iterated maps (among them the 3x + 1 map 
and the related Collatz conjecture, see below) Likely, he should been amused 
to know that there exist simple CA which "compute" this map. Unfortunately, 
only a very faint light seems to arise from these CA concerning the proof of the 
Collatz conjecture. 

2 The 3x+l map and the related conjecture 

The 3x + 1 iterated map was introduced in 1937 by L. Collatz, it was investi- 
gated by a lot of people and it is also known as Collatz map (sequence), Hasse's 
algorithm, Syracuse algorithm; the related conjecture is known as 3a: + 1 prob- 
lem, Collatz conjecture, Ulam problem, Kakutani's problem, Syracuse problem, 
Thwaites conjecture (a large literature on this map, on the related conjecture 
and on possible generalizations is available, see [3-16]). 
The Collatz map (CM): 

u(t + 1) = 3tt(f) + 1 if u(t) is odd; (la) 

u(t + 1) = — ^— if u(t) is even; (lb) 

t € N, u(0) positive integer . (lc) 
The Collatz conjecture (CC): 

For any initial u{Q),the CM attains the final cycle ...1,4,2,1 ... (2) 

Surprisingly enough, the number of steps required to attain the final cycle 
varies sensibly (and apparently unpredictably) also for contiguous small initial 
numbers: f.i. 

u(0) = 6^ 3,10,5,16,8,4,2,1...; 

w(0) = 7 -> 22, 11, 34, 17, 52, 26, 13, 40, 20, 10, 5, 16, 8, 4, 2, 1...; 
u(0) = 8 -> 4,2,1...; 

thus respectively 8, 16, 3 steps are needed. Of course numerical computations 
of CM always confirmed CC, but, in spite of the simplicity of the CM itself, no 
proof (or disproof) of the CC is yet available. Lagarias (1985) showed that there 
are no nontrivial cycles with length <275000 [3,4]. Conway (1972) proved that 
Collatz-type problems can be formally undecidable [5]. Thwaites (1996) has 
offered a £1000 reward for resolving the conjecture [Tl). 
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3 Two CA that "compute" the Collatz map 

Let us now introduce two simple Cellular Automata that "mimic" the CM. 
Looking at the CM itself it is clear that computations should be very simple 
in basis 2 or 3: we take advantage of this. Indeed the first automaton (CA2) 
computes the map in basis 2, while the second one (CA3) computes it in basis 
3. 



3.1 CA2 

This is an unidimensional CA, the cells are arranged on a line, they can be 
numbered by n 6 Z and each cell at the discrete time t (t = 0, 1, 2, ...) can be 
in one of two different states so that the state-function u(n, t) takes values in a 
finite set: u(n,t) € Z/2Z , say {0, 1} .The vacuum-state is 0, i.e. 

u(n,t) — ► as \n\ — ► oo. (3) 

Example: 

...00011010000111010001011001000... 

The evolution law, that allows to construct the state u(n, t + 1) of the CA 
from the known u(n,t), consists in the following rules: 

rule Ia: running from the left to the right (increasing n) tag the cells from 
the first followed by a 1 to the first followed by (extrema included); 

applying this rule in the above example, and tagging with a superimposed 
dot, we get: 

...oooiioioooomoioooioiiooiooo... 

rule Ib: tag also, if the case, from the next 1 following a 1 to the first 
followed by (extrema included); 

applying this rule in the above example we get: 

...oooiioioooomoioooioiiooiooo... 

rule Ic: proceeding to the right, apply rule Ib whenever possible: 
applying this rule in the above example we get: 

...oooiioioooomoioooioiiooiooo... 

rule II: the evolved state obtains according to: 
u(n—l,t+l) = (u(n,t)+u{n-\-l,t)) mod 2. if the cell n is not tagged; (4a) 



u(n— 1, t+1) = (u(n, t) + u(n + 1, i) + 1) mod 2, if the cell n is tagged. (4b) 
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Applying this rule in the above example we get: 
...00100010010100010111001110000... 



Let us now show how this CA computes the Collatz map. Given a CA2 
configuration u(n, t), f.L: 

u(n,t) = ...00010011000..., 

consider the configuration u(n, t) which goes from the leftmost 1 to the 
rightmost 1: 

u{n,t) = 10011. 

A configuration u(n,t) is initially finite, due to ©, and stays finite during 
the evolution, due to rules I and II. Thus we can interpret u{n, t) as an integer 
positive (odd) number in basis 2 (proceeding from the left to the right): 
u(n, t) = 10011 = 1 * 2° + * 2 1 + * 2 2 + 1 * 2 3 + 1 * 2 4 = 25 
It is straightforward to see that, considering the configurations u(n,t) as 
positive (odd) integers, CA2 emulates l(Ta)l. 

Proof. To multiply a number by 3 (namely 11 in basis 2), we have to perform 
the binary summation of that number and of the same number times 2 (which 
in our CA2 configuration is the number itself shifted to the right). In other 
words we have to add any digit to the following one, taking care of the possible 
amount to be carried: indeed we have to carry 1 after a 1 following a 1 up to 
the first following a 0. This justifies Rules lb, Ic, and II. Rule la just adds 1 
to the result. ■ 



Now note that the CM rule l)lb|) in our framework is just a trivial shift of 
the configuration. Thus we can assert that this CA emulate the Collatz map. 
The Collatz conjecture (0) in the language of this CA reads: 

"any initial (non zero) configuration evolves attaining eventually the stable 
configuration ...0001000..." 

Let us give a simple example: 
u(n,Q) = ...00010011000... 
u(n,l) = ...00011001000... 
u(n,2) = ...00101110000... 
u(n,3) = ...00011010000... 
u(n,4) = ...00100010000... 
u(n,5) = ...00101100000... 
u(n,6) = ...00010100000... 
u(n,7) = ...00000100000... 

Interpreting now the configurations as integers, we have the sequence: 
25,19,29,11,17,13,5,1... 

which is the odd subsequence of the Collatz sequence: 

25,76,38,19,58,29,88,44,22,11,34,17,52,13,40,20,10,5,16,8,4,2,1... 

Note the efficiency of CA2 which reaches the final states in only 7 steps 

(versus the 23 steps needed by the Collatz map). The typical behavior of this 

automaton is shown in Figure 1. 
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3.2 CA3 



This is an unidimcnsional CA, the cells are arranged on a line, they can be 
numbered by n £ Z and each cell at the discrete time t (t = 0, 1, 2, ...) can be 
in one of four different states so that the state-function u(n, t) takes values in a 
finite set: u(n, t) £ Z/2Z , say {0, 1, 2, 3} .Moreover 

u(n,t) — > ; (5a) 

n — > — oo 

u(n,t) — ► 3; (5b) 

n — > + oo 

If u(n,t) = 3 ift,en u(n + l,t) = 3. (5c) 
...00012010222102110012333... 

The evolution law, that allows to construct the state u(n, t + 1) of the CA 
from the known u(n, f), consists in the following two rules: 

RULE III: tag the cells following a non tagged 1 up to the next 1 or 3 
(included); 

applying this rule in the above example, and tagging with a superimposed 
dot, we get: 

...0001201022216210110012333... 

rule IV: the evolved state of CA3 obtains according to the following table: 



tagged cells not tagged cells 

-> 1 0^0 

1 -»• 2 1 -> (6) 
2^2 2 -> 1 

3^2 3 -> 3 

applying this rule in the above example we get: 
...0000212011101220020002233... 

Let us now show how this CA computes the Collatz map. Given a u(n,t), 



f. 



u(n,t) = ...000101121333..., 

consider the configuration u(n, t) which goes from the leftmost not-zero num- 
ber (included) to the first 3 (excluded): 
u(n,t) = 101121. 

A configuration u(n,t) is initially finite, due to (JSJ, and stays finite during 
the evolution, due to rules III and IV. Thus we can interpret u(n, t) as an integer 
positive number in basis 3 (proceeding from the right to the left): 

u{n,t)= 101121 (basis 3)= l*3° + 2*3 1 + l*3 2 + l*3 3 + 0*3 4 + l*3 5 = 286 
(basis 10). 
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Suppose now that u(n, t) is even (as in the above case); then in u(n, t) there 
is an even number of 1 and therefore, according to Rule III, in u(n, t) the first 
3 cannot be tagged: 

u{n,t) = ...000101121333..., 

In this case Rules III and IV perform just the division of the number u{n, t) 
by 2 , according to the CM rule Ijlbjl. The proof of this statement is straight 
forward and it is based mainly on the trivial fact that 

3 k+1 = (2 + 1) *3 fe . (7) 

Let us explicit this using the above example: 

u(n, t)= 101121 = 1 * 3 5 + 0*3 4 + l*3 3 + l*3 2 + 2*3 1 + l*3° = 

(2 + 1) * 3 4 + * 3 4 + 1 * 3 3 + 1 * 3 2 + 2 * 3 1 + 1 * 3° = 

2 * 3 4 + (2 + 1) * 3 3 + 1 * 3 3 + 1 * 3 2 + 2 * 3 1 + 1 * 3° = 

2 * 3 4 + 4 * 3 3 + (2 + 1) * 3 1 + 2 * 3 1 + 1 * 3° = 

2 * 3 4 + 4 * 3 3 + * 3 2 + 4 * 3 1 + (2 + 1) * 3° + 1 * 3° = 

2 * 3 4 + 4 * 3 3 + * 3 2 + 4 * 3 1 + 4 * 3°. 

Thus clearly: 

u(n,t)/2 = (basis 3)=12022 (basis 3)= 143 (basis 10). 

Indeed the tagging Rule III and the Rule IV perform this task on the CA 
configuration : 

u(n,t) = ...000101121333...; 
apply Rule III: 
...000101121333..., 
apply Rule IV: 

u(n, t + l)= ...000012022333... ii(n, t + 1) = 12022. 
Let us now consider the case of odd u(n, t), f.i. 

u(n,t) = ...000101120333... u(n,t) = 101120 (basis 3)=285 (basis 10). 
CM rule JTaJ yields 

285 * 3 + 1 (basis 10)=856 (basis 10)=1011201 (basis 3). 

It is plain that in terms of the CA3 configuration this could be accomplished 
just by changing the leftmost 3 into a 1. But then the new number is even 
and thus one should again apply CM rule (|lb(l yielding as a result: 

285 => 856 => 428 (basis 10)=120212 (basis 3). 

It is straightforward to see that Rules III and IV perform these two consec- 
utive tasks in just one step. 
In the above example: 

...000101120333... Ru J^4 u ...000101120333... Ru M4 v ...000012021233... 

Note that, even if all "computations" are in basis 3, we need four states for 
CA3. Another oddity, in terms of usual unidimcnsional CA, is the presence of 
two different vacuum states (0 on the left, 3 on the right). 

Note also that the final cycle, eventually attained according to the Collatz 
conjecture, is visualized by CA3 in the following way: 

...0001333... Ru M4 n ...0001333... Ru M4 v ...0000233... Ru J^4 n ...0000233... 1 
...0000133... 

A typical behavior of this CA is shown in Figure 2. 
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4 Final remarks 



Looking at Figures 1 and 2, we can see that one edge of the configuration moves 
with constant average velocity while the other exhibit erratic behavior: this 
somehow justify the conjecture ( obviously prove not!). However if the intro- 
duction of these two CA doesn't seem to give a great contribution to the aim of 
proving the Collatz conjecture, on the other hand it allows to explore intensively 
the 3x + 1 map also for very large initial integers. Thus our CA offer a good 
example of how CA themselves can be used as fast computing machines even 
if implemented on sequential computers (much more if running on dedicated 
parallel machines). Moreover this approach allows convenient investigation of 
the statistical behavior of iterates. Indeed in J7|, the authors proposed a Quasi 
Cellular Automaton for the 3a; + 1 map: they did not introduce a real CA but 
merely displayed the iterated of CM in basis 2, simulating in fact our CA2. 
In the same work generic behaviors associated with classes of seed values and 
periodic and chaotic structures within iterate patterns were investigated. The 
same investigation could be pursued for our CA3. 
Finally some further remarks: 

• our CA, according to the Wolfram heuristic classification JH], should be 
Class 1 CA (thus they should be not interesting ones...); 

• our CA, due to the tagging rules I and III, are "non local" CA (namely 
the evolved state of a cell depends in principle on the whole previous con- 
figuration) . Indeed the evolution laws for these C A exhibit some analogies 
with the so called "fast rule" introduced in for the parity-rule "filter" 
CA [20| : this suggests that it should be possible (may be useful) to find a 
" local" , possibly " filter" , formulation of these C A; 

• empirically, computer investigation shows that "small deformations" of 
the evolution law of CA2 and, to a lesser extent, of CA3 give rise mostly 
to new CA which themselves end up with final cycles (see Figures 3,4,5: 
of course this is not equivalent to change " slightly" the rules of the Collatz 
map, thus we have non conventional generalizations of the map itself); 

• the scheme here introduced could be easily and profitably extended to 
emulate and explore other different iterated maps. 



Figure captions 

Figures 1-5 are space-time diagrams of the evolution of the involved CA, 
space on the vertical axis, time flows left-right on the horizontal axis, the states 
of the cells are shown in different shades of grey. 

• Figure 1: the typical behavior of CA2, starting from a random chosen ini- 
tial configuration of 300 cells (corresponding to an integer of order 2 "300) 
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and ending with the final cycle predicted by the Collatz conjecture, is 
shown. 

• Figure 2: the typical behavior of CA3, starting from a random chosen ini- 
tial configuration of 100 cells (corresponding to an integer of order 3" 100) 
and ending with the final cycle predicted by the Collatz conjecture, is 
shown. Note that, in order to get a final stable configuration, the average 
drift velocity was corrected, adding a downward shift at any 2 time-steps 

• Figure 3: the typical behavior of a variant of CA2, starting from a random 
chosen initial configuration of 450 cells (corresponding to an integer of 
order 2~450) and ending with the same final cycle of CA2, is shown. The 
evolution law of this CA is the same of CA2 except for the Rule lb that 
now reads: tag also, if the case, from the next 1 following a 1 or followed 
by a 1 to the first followed by (extrema included). 

• Figure 4: the typical behavior of a variant of CA2, starting from a random 
chosen initial configuration of 400 cells (corresponding to an integer of 
order 2 '400) and ending with a final cycle (more complex than that of 
CA2), is shown. The evolution law of this CA is the same of CA2 except 
for the ending condition of tagging in Rule la, lb that now reads: ... to the 
first followed by and following a (extrema included). 

• Figure 5: the typical behavior of a variant of CA3, starting from a random 
chosen initial configuration of 200 cells (corresponding to an integer of 
order 3~200) and ending with a simple final cycle, is shown. The evolution 
law of this CA is the same of CA3 except for Rule IV, that now reads: 
tag the cells following a non tagged 1 up to the cell that is 1 or 3 or is 
followed by a 2. 
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